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Abstract 

We apply the Kovacs experimental protocol to classical and quantum p- 
spin models. We show that these models have memory effects as those ob- 
served experimentally in super-cooled polymer melts. We discuss our results 
in connection to other classical models that capture memory effects. We pro- 
pose that a similar protocol applied to quantum glassy systems might be useful 
to understand their dynamics. 
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1 Introduction 



The Kovacs memory effect was first reported by this author in the 60s pQ. It 
demonstrates that super-cooled liquids when taken out of equilibrium have a very 
intricate dynamics that cannot be predicted on the basis of the knowledge of the 
instantaneous value of the state variables (P, V, T) right after the perturbation. More 
precisely, Kovacs showed that the specific volume of a polymer melt in its super- 
cooled liquid phase, has a rather non-trivial evolution that depends on the thermal 
history of the sample PQI21- Non-trivial effects of temperature variations were also 
observed in the evolution of (two-time) susceptibilities of dipolar glasses spin- 
glasses |lj and many other glassy systems El- History-dependent phenomena 
in granular compaction have been recently reported [7, . In this case, the control 
parameter is the tapping strength and the observable is the density. 

The experimental setup involved in the Kovacs effect is the following. In a first 
step, one quenches an equilibrated liquid with a very fast cooling rate from a high 
temperature T to a low temperature T 2 , at time t — 0. One then follows the 
subsequent evolution of the quantity of interest that in Kovacs' experiments is the 
specific volume. In order to match what we shall study in this paper, we describe 
the experiment using another one-time quantity, the energy density, as the example. 
The energy density relaxes in time and it slowly approaches an asymptotic value 
that may fall out of the experimental time-window. The time-dependent curve 
£( T2 \t) constitutes a reference and we use a superscript (T 2 ) to indicate that it has 
been obtained using the first prescription, at fixed external temperature T 2 . The 
extrapolated asymptotic value defines £ as (T 2 ). In a second step, one first quenches 
the sample from the same high temperature T to a lower temperature Ti(< T 2 ) 
at time t — 0, waits until a carefully chosen time ti and then heats the sample 
to the final temperature T 2 . The value of t\ is chosen such that £ <yTl ^ T2) {ti) = 
£ as (T 2 ), where the plus sign indicates that the one-time quantity should match the 
asymptotic value obtained with the first procedure right after heating the sample 
from Ti to T 2 . [In the experimental protocol the need to use £^ Tl ^ T2 \tf) = £ as {T 2 ) 
is due to the fact that when changing the temperature there is a trivial response 
due to the thermal expansion of the local degrees of freedom.] Since the initial 
value of the energy density at the final temperature T 2 is already the asymptotic 
one, one could have expected that the energy remained constantly fixed to this 
value for all subsequent times independently of the value of T\ (apart from very fast 
rearrangements decaying exponentially). However, Kovacs demonstrated that after 
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ti the energy-density, £^ Tl ^ T2 \t) , has a slow non-monotic dependence on time, first 
increasing and then decreasing back to its initial and asymptotic value £ as (T 2 ): 

g( T ^\t) = £ as (T 2 ) + A£(t) , (1) 

with the "Kovacs hump" A£ satisfying 

A£(t) > , A£(tf) = lim A£(t) = . (2) 

t — >oo 

The form of the hump depends on the values T 2 and T% used. Qualitatively, its 
height is an increasing function of T 2 — Ti and the time at which the maximum is 
reached decreases when T 2 — T% increases. 

It is interesting to note that Kovacs' experiments have been performed in the 
super-cooled liquid phase. While it would be no surprise to find nonequilibrium 
effects in the glassy phase, the reason why one also finds a nonequilibrium behaviour 
in the super-cooled liquid is that the jump in the external temperature drives the 
system out of equilibrium and the relaxation occurs in a very long time-scale. This 
experiment proves that the knowledge of the state variables, in Kovacs' experiments 
P, V and T at tf, is not sufficient to determine the subsequent evolution of the 
same quantities in the glassy and super-cooled liquid phases (if the latter has been 
recently strongly perturbed). 

Recently, several authors presented analytical and numerical studies of this effect 
using a variety of models with glassy dynamics. So far, apart from phenomenological 
approaches |2J, the Kovacs effect has been analysed numerically with molecular 
dynamic simulations of a molecular model of a fragile glass former |3 and Montecarlo 
simulations of the 3d Edwards- Anderson spin-glass [5], and analytically within the 
ferromagnetic Ising chain [TO] , the critical 2d XY model the trap model |12j . 
domain growth [12], 1^ kinetically constrained spin models of fragile and strong 
type [13], and the parking lot model of granular matter [14 . It has been suggested 
that the quantitative analysis of the Kovacs effect may help distinguishing between 
different glassy models and, perhaps, may also help identifying spatial properties 
of glassy systems. In this paper we show that the main qualitative features of the 
classical Kovacs effect are captured by the fully connected spherical p-spin disordered 
system [15] . a model with no spatial structure but with a slow dynamics leading to 
very slow relaxations close to the transition to and in the glassy phase JHl E]- 
This 'negative' result, as far as what can be deduced about spatial rearrangements 
from the Kovacs effect, is reminiscent of the discussion ^H] on the interpretation of 
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hole burning experiments JH]- We also discuss the scaling laws that describe the 
behaviour of the hump and compare them to what found in other glassy models. 

On the other hand, the analysis of quantum glassy systems is now starting to 
call the attention of experimentalists and theoreticians. The slow, history-dependent 
relaxation of a dipolar quantum system in its glassy phase has been reported [20J. 
The sample that entered the glassy phase following a quantum route (changing the 
strength of quantum fluctuations at fixed low temperature) is always in advance 
with respect to the one that arrived at the same point in parameter space following 
a classical path (keeping the strength of quantum fluctuations fixed and cooling the 
system) . This has been demonstrated by the fact that the time- dependent dielectric 
constant of the quantum-cooled sample is closer to its asymptotic value at any 
finite, experimental time. Memory effects were also observed in glasses at ultralow 
temperatures |2I] and the electron glass 

A variant of the Kovacs procedure where the control parameter is the strength 
of quantum fluctuations can be easily envisaged. The question then arises as to 
whether a hump appears and which is its structure, scaling form, etc. We address 
this question using a quantum extension of the p-spin model introduced and studied 
in [231E1! (see also |25|) - 

In short, in this paper we show that simple disordered mean-field models capture 
the phenomenology of the Kovacs effect. With this aim, we analyse the nonequi- 
librium relaxation of the spherical p-spin disordered model in its classical [To^ IT7] 
and quantum versions, see Sect. |2 for their definitions. First, we reproduce the 
classical setup using temperature as the control parameter and we discuss the re- 
sults in comparison with previous explanations of the same effect (Sect. EH). Second, 
we switch on quantum fluctuations and use their strength as the control parameter 
(Sect. 0J). In both cases we follow the evolution in time of the potential energy- 
density of the system. We present our conclusions in Sect. E3 

2 The spherical p-spin model 

The spherical p spin model is defined by the Hamiltonian [T5] 

Hj[S] = — ^ Jiii 2 ...i p s ii s i2 ■ ■ ■ s i p ■ (3) 
(ili2...i p ) 

The spins s« are continuous variables Sj G (—00,00) forced to satisfy the global 
spherical constraint Yli=i s i — N with N their total number in the sample. The 
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exchange interaction Ji x i 2 ...i are quenched random variables drawn from a Gaus- 
sian distribution with average [/^...ij = and variance [J^...^] = J 2 p\/(2N P ^ 1 ). 
We henceforth use square brackets to indicate the average over disorder. The in- 
teractions occur between all groups of p spins in the sample. The model is then 
fully- connected and mean-field in character. The scaling of the variance [J? i ] with 
N has been chosen so as to ensure a good thermodynamic limit. The parameter p 
takes integer values, p > 2. Mixed models with a Hamiltonian with two terms of 
the form @ with p = 2 and p = 4 |2Ej are also of interest [Tj\ . 
The classical dynamics is determined by the Langevin equation 

Si(t) = -^^-»(t)si(t)+Zi(t) (4) 

where E(t) is the total energy, /i(t) is the Lagrange multiplier enforcing the spherical 
constraint, and is a Gaussian white- noise with zero mean and correlations 

m)Ut')) = ^ B T5 l3 5(t-t') . (5) 

We assume that the system is prepared at t = with an infinitely fast quench from 
T — > oo to the initial temperature T\. The dynamic equations for the macroscopic 
order parameters are derived using standard functional techniques [T7]. We discuss 
them below, as the classical limit of the quantum extension of the same model. The 
models with p = 2, p > 3 have different dynamics; the former yields a mean-field 
description of simple domain growth, while the latter case is related to the schematic 
mode-coupling theory (mct) of super-cooled liquids and glasses. These models and 
their generalizations are reviewed in detail in [TTj . 

Quantum fluctuations can be introduced [2S] by upgrading the spins to coordinate- 
like operators, adding a 'kinetic term' 

1 N - 

— Ytf (6) 
2M f-i 1 K J 

to the Hamiltonian 0, and imposing canonical commutation relations 

[Ilj, sj] = ih 5ij , [s^ sj] = , [Ilj, Uj] = . (7) 

At time t = we set the model in contact with an Ohmic bath of quantum harmonic 
oscillators 

N i N i 

H b = Ett— P? + Eo m ^?. (8) 

^6 = - E «i E c ^ £ ^ > ( 9 ) 

i=l i=l 
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where, for simplicity, we considered a bilinear coupling, H^. (Note that for the 
spherical problem it is not necessary to introduce a counterterm). We assume that 
this environment has a well-defined temperature and that it is not modified by 
the interaction with the system. The initial density matrix is factorised and we 
choose random initial conditions for the system. After integrating out the bath 
degrees of freedom, and using the fully-connected character of (j3J), one arrives at a 
dynamic generating functional from which one derives exact dynamic equations for 
the macroscopic two-time dependent order parameters 



i^r/r.m ft ,„vm D ,wx 1 £ [<*(*)>] 



■ (io) 

h=0 



is an infinitesimal field that couples linearly to the spins, modifying the Hamilto- 
nian as H — > H — J2i=i hi§i. C is the symmetrized correlation function and R is the 
linear response of the system. 

The dynamic equations take the Schwinger-Dyson form 

[Md* + fi(t)}R(t,t') =8{t-t')+ f dt"m,t")R{t",t') , (11) 

J a 

[Md'? + ii(t)]C(t,t')= /* dt"m,t")C(t",t')+ f dt"D(t,t")R(t',t") ,(12) 

Jo Jo 

with the self-energy S and the vertex D given by 

£(*,*") = -A-qit-t') +a(t,t') , (13) 
D(t,t") = 2hu{t-t') + d(t,t') . (14) 

The first contributions originate in the interaction with the Ohmic bath of spectral 
density [23121! 

I(w) = —oje^eiu) , (15) 

7T 

where A is a cut-off included to avoid divergences. The kernels rj and u are functions 
of the time-difference r = t — t' and they are given by 

4 7 A 2 2Ar , . 

viT) = (l + (Ar)T' (16) 



-- ^ I" du cue-^coth f ^ ) cos(cur) . (17) 

7T JO \ 2 J 

The second contributions, a and d, are due to the interactions in the system and 
read 



cr(t,t') = — Im 



vJ 2 ih 1P ~ 1 

>■' " C(t,t')-^R(t,t') 
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d(t,f) 



pJ 2 



Re 



C(t,t')-^(R(t,t')+R(t',t)) 



p-i 



(19) 



An integral equation that fixes the Lagrange multiplier fi(t) supplements these equa- 
tions and it is derived from the requirement C(t,t) = 1 for all times [2*3] . 

The classical limit is easily obtained by neglecting the kinetic term and by taking 
the limit h — > 0. Indeed, in this case, the effect of the coupling to the bath reduces 
to the usual contributions originating in the friction and noise terms of the Langevin 
equation and involving a first-order time-derivative of the correlation and response 
when A is further taken to infinity. 

In the quantum case, three contributions to the total energy density can be 
identified: the kinetic part, the potential part and the interaction with the bath. In 
the following we focus on the averaged potential energy density 



E(t) = N-*[(Hj[S\)] 



1 

'N 



^ Jili2—i p Sii (0**2 (0 • • • Si p {t) 
v(iii2...* p ) ' 



(20) 



With a simple calculation one finds 

roo 

Bit) = / dt' [a(t, t')C(t, t') + d(t, t')R(t, t')} 
Jo 

that in the classical case becomes 

T 2 r> rt 



(21) 



S(t) = -— [ dt' C p -\t,t')R{t,t') . 
2 Jo 



(22) 



3 Classical case 

In this Section we analyze the classical problem with p = 3. After recalling the value 
of the asymptotic energy-density for the isothermal relaxation, we solve numerically 
the dynamic equations for C and R following Kovacs' protocol. We analyse the 
Kovacs hump and we obtain its scaling with time and temperature in the three 
regimes of short, intermediate and long times. 

3.1 Analytic results 

In the high-temperature phase the Fluctuation Dissipation Theorem (fdt) implies 
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ks = J = 1 henceforth. We shall focus on the p = 3 model and use T 2 = 0.75 as 
the classical reference case for which £ as (T 2 ) ~ —0.67 (see Fig. [TJ. 

This model undergoes a dynamic transition from an equilibrium (paramagnetic) 
to a nonequilibrium (glassy) phase at 



P(P - 2 ) P - 2 (24) 



For p = 3, Td ph 0.61 and the asymptotic energy- density at the dynamic critical 
temperature takes the value £ as (Td) ~ —0.82. 

In the low-temperature phase the solution of equations (fTTj) -(lT9 j) involves a modi- 
fication of the FDT [iniQZI, Rs(t,t') = T~/ f d tl C s (t,t')9(t-t'), for the slow part of the 
relaxation. The effective temperature is given by T e ff = (q ea T)/[(p — 2)(1 — q ea )} 
and the energy- density reads 

£as{T) = ~ [(p - 2)(1 - q ea )q£ 1 + (1 - q p j] (25) 
with the Edwards-Anderson order parameter, q ea , determined by 

^^e 2 (l - leaf = T 2 . (26) 

At the dynamic transition T e ff = T, q^ = (p — 2)/(p— 1) and is given by Eq. pijl . 
The numerical solution of Eqs. (|2*Hj) and (j2B| yields the value of the asymptotic 
energy- density in the glassy phase. 

3.2 Numerical results 

The effect of temperature variations on the dynamics of the p spin models in the 
glassy phase was studied in a couple of papers. The effect of small amplitude tem- 
perature cycles on the nonequilibrium relaxation of the p = 2 model were discussed 
in |2H] while their influence on the dynamics of the p = 3 was analyzed in • 

We show here that the p spin model with p > 2 captures a similar phenomenology 
in the sense that a hump with a slow relaxation is obtained when the Kovacs' protocol 
is applied. The model with p = 2 is related to ferromagnetic domain growth, as 
described by the O(N) model in d — 3, and our results are intimately related to 
the ones in ^21 E] • In the following we present the data for p = 3 (the schematic 
mode-coupling-theory [HH^H]) onr y- 

Figure Q shows the time-evolution of the energy-density in the p = 3 classical 
model using different temperature jumps chosen according to Kovacs' rule. The 
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Figure 1: The evolution in time of the energy-density £(t) in the classical p = 3 
spin model close and above the paramagnetic - spin-glass transition. The solid line 
is £^ T2 \t) and has been obtained using a rapid quench to T2 = 0.75. The other 
curves show £( Tl ~* T3 ' (t) , i.e. the result of having performed a temperature jump 
from Ti = 0.65, 0.7, 0.725 to T 2 at ti = 92, 170, 295 time-steps, respectively. The 
asymptotic value £ a s{T2) ~ —0.677 is shown with a horizontal line. 

solid line has been obtained using a quench to T 2 = 0.75. Note that since T 2 > Td 
the quench is done within the paramagnetic phase in which the equilibration time 
is finite. The asymptotic value approached with the algorithm is S as (T 2 ) ~ —0.677. 
(The analytic prediction is £ a s(T 2 ) ~ —0.667. Here and in what follows we use 
a time step 5 = 0.05 to numerically solve the integro-differential equations.) The 
other curves include a temperature jump from Ti to T 2 at t\ when S^ Tl ^ T2 \ti) = 
^asiT-i). The values of the parameters Ti and t\ are given in the caption. As 
observed by Kovacs experimentally PP and for a variety of models jH] [E|, the 
evolution of the energy-density is non-monotonic. Since the time t\ is chosen so 
as to have £ Tl ^ T2 (t^) = £ as (T 2 ) ~ —0.677, after the temperature jump the energy 
density already has the expected asymptotic value. However, it evolves during 
approximately three decades before reaching the reference curve, first increasing 
towards a maximum and then decreasing towards S as (T 2 ). 

The qualitative features of this non-monotonic behaviour are better described 
by analysing the hump AS. Figure 121 shows its dependence with the temperature 
difference T 2 — Ti. As in the experimental data, the height of the hump, ASk = 
max AS, is an increasing function of x = T 2 — Ti, while the position of the maximum, 
tx = t\ max , is a decreasing function of the same parameter. In Fig. El we show that 
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Figure 2: The hump, A£(t), as a funtion of the time-difference, t — t\ 1 in the classical 
p spin model with p = 3 close and above its paramagnetic - spin-glass transition. 
The different curves correspond to the intermediate temperatures T\ = 0.4, 0.5, 0.6, 
0.65, 0.7, 0.725 (from top to bottom) and the reference curve has been obtained for 
T 2 = 0.75. 



these quantities are well described by functions of the type 



AS K (x) 
t K (x) 



ax (1 — b In x) 
— (1 — dlnx) 



(27) 
(28) 



x 



with a, 6, c, d fitting parameters. 

We also show in Fig. 01 that tx, the position of the maximum, can be fitted as 
tx = oi\, where t\ is the time when the temperature jump has been applied, and a 
and b are constants. 

Three time-regimes can be identified in the hump: short times well before the 
maximum is reached, intermediate times around the maximum and long-times when 
the hump approaches its asymptotic vanishing value. The time-temperature depen- 
dence in the three regimes can be summarized as follows. 

Short times 

In the short time regimen, a linear function with logarithmic corrections ^2] 



fits the data with great precition. Indeed, in Fig. El we show together the data for 
T\ = 0.6, 0.65, 0.7 and t\ = 63, 92, 170 respectively and the fits by J- S (x). 



T s {x) = a x (1 — b In x) 



(29) 
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Figure 3: The maximum height of the hump, A£ic(t) — max A£ (red curve, +) and 
its position tx = t\max (green curve, x) as a funtion of the difference in temperature 
T2 — T\ in the classical p spin model with p = 3 at T2 = 0.75. The lines are fits 
to power laws with logarithmic corrections, see Eqs. (|27|l and J5SJ), with parameters 
a = 0.153 and b = 0.556 for the height of the maximum and c = 16.73 and d — —0.135 
for its position. 
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Figure 4: The position of the maximum in the hump, Ir = t\ max , as a function 



of the time when the jump in temperature is imposed, t%. 
relation t K = at\ with a = 1.901 and b = 0.908. 



The line is the power-law 
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The data can also be scaled using 



A£(t) = (r 2 -ri)^.((t-ti)^) 



(30) 



An accurate description of the rescaled data is achieved taking a = 0.0623 and 
b = 0.251 for F s {x), see Fig. El 

Intermediate times 

The time-temperature dependence of the intermediate part of the relaxation is 
well described with the scaling form 



see Fig. |7J This scaling law is of the class found in EH • 
Long times 

The long times decay of the hump, i.e. for times well beyond fx, is exponential 



The curves AS(t) for different T\ can be made to collapse by shifting time according 
to 



that is to say that, the curves for the systems on which a temperature shift was 
applied are in advance with respect to the reference one. 

Finally, we implemented the same protocol using temperatures T\ and Ti that are 
both below the dynamic critical temperature As shown in Fig. [01 the qualitative 
behaviour is the same as on the other side of the transition point. It is interesting 
to note that the curves for the perturbed system join, and become independent of 
T\ before reaching the reference curve (this is similar to what observed in (Hj). 

In the low temperature phase the scalings discussed above are modified as follows. 
In the short-times regime, the scaling function ()29)1 remains valid but the scaling 
form is modified to 




(31) 



A£(t) = ae 



-bt 



(32) 



(33) 




(34) 



a 



T 2 -T : 



(35) 
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Figure 5: The short time fitting of the data for three values of the lower temperature, 
Ti = 0.6,0.65,0.7, using Eq. <P|). 
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Figure 6: The short time scaling of the data for three values of the lower tempera- 
ture, Ti = 0.6,0.65,0.7, using Eq. (J3J3J. The solid line is a fit of the rescaled data for 
Ti = 0.65 using Eq. with a = 0.0623 and b = 0.251. 
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Figure 7: The scaling for intermediate times with Ti = 0.5, 0.6, 0.65, 0.7, see the 
text for an explanation. 
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Figure 8: For long times the curves collapse on a master curve under the time-shift 
in Eq. i|33|) . With red solid line the reference curve £^ T2 \t); with dashed lines the 
curves under the effect of the temperature jump from T\ = 0.5, 0.6, 0.65, 0.7 translated 
in time t — > t + At with At = 15, 19, 23, 31. The final approach to zero is exponential. 
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Figure 9: The hump in the glassy phase of the classical p-spin model with p = 3 
using Ti = 0.3, 0.4, 0.45, 0.5 with t x = 80, 173, 330, 1080, respectively, and T 2 = 0.55. 



For intermediate times we found 

A£(t) = A^((A^-) 1+Q ) , (36) 

with Ti independent of Ti and T 2 . At very long times, the approach to zero is a 
power law, at~ b . The shift t — > t + (tx — ti) is also efficient collapsing the data. 



3.3 Discussion 

The non-trivial content of the hump is its slow relaxation and dependence on T<i—T\. 
The observed behaviour has been captured by several models already presented in 
the literature. Let us discuss these explanations and compare them to the one 
associated to the p spin model. 

First, one can compaire the behaviour observed in the super-cooled liquid phase 
with a simpler system: an overdamped harmonic oscillator. For this model a direct 
calculation yields 



AS(t) = - 
v ; 2 



k(x 2 ) - k B T x ( 1 - ef *i ) - k B T 2 e^' 



(37) 



with k the harmonic constant, 7 the friction coeficient and Xq the initial condition. 
The time t\ is fixed by requiring that A£(ti) = 0. This condition forces the bracket 
in equation (|37j) to be zero and then, there is not a forthcoming hump. Clearly one 
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needs to go beyond this model to get the observed two-temperature dependence and 
slow relaxation. 

A simple next step is to study a model with a distribution of relaxation times that 
depends on temperature. A simple realization is the 2d XY model in the spin-wave 
approximation taking into account the contribution of vortices [TT]. This model 
is given by a Gaussian free field (the angle of the local magnetization) with a In- 
dependent stiffness, p(T). It has a slow dynamics characterised by the growth of a In- 
dependent correlation length ir{t) = (p(T)t) 1 / 2 and it captures the phenomenology 
of the Kovacs' effect, as discussed by Berthier and Holdsworth [TT] . 

In slightly more general terms [3 E] the Kovacs' effect can be rationalized 
in any system with a growing temperature-dependent dynamic correlation length, 
that is shorter than the equilibrium one. When one shifts the temperature to T 2 
at t\, the length scales that are shorter than #ft(£i)> and are hence equilibrated 
at Ti, have to reequilibrate at T 2 , where their equilibrium energy is higher. The 
structure reached at Ti has to break up and allow for the nucleation of new structures 
equilibrated at T 2 . Instead, the length scales that are longer than ^Ti(^i) are still 
not equilibrated at Ti and they may continue their evolution to equilibrate now 
at T 2 . The former processes involve shorter length-scales and should be faster and 
dominate the first part of the relaxation after t\ hence leading to an energy increase. 
The latter processes are slower and dominate the decay from the maximum towards 
the asymptotic value S as (T 2 ). Within this picture, the time at which the maximum in 
AS is reached corresponds to the time when the small length scales have equilibrated 
at the new temperature T 2 . A similar argument was put forward to explain the 
overshoot observed in the time-dependent dielectric constant of dipolar glasses after 
a temperature jump It is also behind the calculation presented by Brawer on 
the ferromagnetic Ising chain at very low temperature fTUl . 

However, it is not necessary to invoke a growing correlation length to capture 
the qualitative features of the Kovacs' effect. The main ingredient in the p-spin 
model that leads to this effect is the slow - non-exponential - and temperature 
dependent relaxation of the linear response after a strong perturbation. A scenario 
with a wide spectrum of relaxation times that depends on temperature was used in 
the past to explain the Kovacs' effect [S]. Here we demonstrated that, as one could 
have expected ^2], the p-spin model or, equivalently, the mct with no equilibration 
assumption, has this property. 

The temperature and time dependence of the hump do depend on the model 
considered but the main qualitative features of the effect are shared by all of them. 
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In the context of the spin models related to the mct these will be obviously modified 
if one considers a mixture of p = 2 and p = 4 models (that corresponds to going 
from the schematic mct to more refined versions). 

Finally, let us discuss the relevance of the Kovacs' effect for the development of 
a thermodynamics of the nonequilibrium glassy state. 

The attempts to use a thermodynamics for nonequilibrium glassy systems are 
based on the introduction of effective state variables, basically temperature and 
pressure. Initially, a (constant) Active temperature that characterizes the glassy 
structure was introduced by Tool [30J. As a consequence of Kovacs' experiments it 
was realized that this single parameter was not enough to describe the evolution of 
the glass and the Active temperature was upgraded to be a full history dependent 
function measuring the departure from equilibrium P^EIE]- Whether the Active 
temperature, as defined in [TJ EJ0J EH] behaves as a thermodynamic temperature 
remains to be proven. One can also translate the study of the Kovacs' effect in 
the parking lot model by Tarjus and Viot ^3] in these terms: the second time- 
dependent state variable introduced generalising Edwards' prescrition [H2] acts as a 
Active temperature. 

More recently, an efiective temperature was deAned using the modiAcation of the 
Auctuation - dissipation theorem in slowly evolving nonequilibrium systems. The 
interpretation of this quantity as a bonafide temperature was discussed, and some 
conditions on the physical relevance of this deAnition were also given [21] . In partic- 
ular, the need to have a system evolving slowly - characterized by a slow relaxation 
of one time quantities - to be able to associate the properties of a temperature 
to the fdt ratio was reckoned and stressed. The situation in Kovacs' experiment 
goes beyond this limit: the system is strongly perturbed and the one-time quantity 
under study has a non-trivial, non-monotonic time-dependence. The Fluctuation 
Dissipation Relation also depends on time in a non-trivial manner jH3]- Thus one 
cannot assert that it leads to a bonafide temperature and use it to construct a simple 
extension of thermodynamics to describe the subsequent behaviour of the system. 

A similar conclusion, though expressed in terms of the potential energy landscape 
(pel) scenario was reached by Mossa and Sciortino [B]. These authors showed, with 
molecular dynamic simulations of the fragile glass former OTP, that two systems in 
identical thermodynamic conditions (same values of T, V, P) can be in very difierent 
regions of their potential energy landscape if one of them has been strongly per- 
turbed. Since the strongly perturbed system wanders in a region of the pel that 
is never sampled in equilibrium, its conAguration cannot be associated to that of 
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an equilibrated liquid at a different temperature. The region of the pel sampled 
allows for a definition of a microcanonical temperature only when the variation of 
the external conditions (temperature in this case) is small. Mossa and Sciortino 
arrived at this conclusion by comparing the properties of the inherent structures 
(closest local minimum in the pel) visited during aging after a temperature jump 
of large magnitude to those sampled in equilibrium. 

4 Quantum model 

When quantum fluctuations are switched on one has to deal with the full Schwinger- 
Dyson equations (|TT|) -(fT^ j) . As explained in (33] the parameter that plays the role 
of the transverse field in trully quantum spin models is here the inverse of the 
mass M. Indeed, this model undergoes a phase transition in the (T, V) plane with 
T = h 2 / (MJ) from an equilibrium paramagnetic phase to a non-equilibrium glassy 
one. In order to test the memory effect in the quantum problem, we then apply the 
Kovacs' protocol using M~ l as the control parameter and we follow the evolution 
of the averaged potential energy density £(t). 

In the quantum problem the asymptotic value of the potential energy density 
depends on M. (In the classical limit it does not.) The relaxation of the potential 
energy density at constant M has (damped) oscillations whose magnitude depend 
on the parameter M. For large values of the mass the oscillations have a suffi- 
ciently large amplitude such that the asymptotic value falls within the oscillation, 
i.e. S(t) < S as (T2,T2) for some finite times t. In the following we choose a value 
of M 2 such that the system is close to the paramagnetic - glass transition and for 
which £ as (T 2 ,r 2 ) < £{t) for all finite time t, see Fig. [TU| Another feature to signal 
is that the oscillation in the energy- density decay may be such that there is more 
than one value of t\ for which £ Ml (tf) = S^ 2 , see the curve drawn with a dashed 
line in Fig. El for which M = 0.8. 

Figure HP shows the result of applying Kovacs' protocol to a system at T = 0.75 
with the reference value of the mass, M 2 = 0.5. We use four values of M 1; Mi = 0.6, 
0.8, 1 and 1.2 that satisfy Mi > M 2 , i.e. Ti < T 2 . For M x = 0.6, 1 and 1.2 we find 
a unique t\ satisfying £ Ml {ti) = £^ 2 . For Mi = 0.8 instead three solutions to this 
equation exist as shown in Fig. Displayed in Fig. ^2 ar e the reference energy 
density £( T2 ' r2 \t) (solid red line) and the hump in the energy densities obtained by 
shifting the mass. 

The first thing to note in the figure is that the curves depend on the value of Mi , 
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Figure 10: Energy-density decay in the quantum model with p = 3 at constant 
mass and temperature. Solid line M = 0.5 and dashed line M = 0.8. In both cases 
T = 0.75. 

similarly to what happened in the classical case with 7\: the larger the difference 
in the masses (or in the quantum parameter T), the more pronounced the effect. 
There is also a weak decreasing dependence of the position of the maximum with 
M 2 - Mi (and T x - T 2 ). 

A second feature to remark is that the curves following the jump in M go above 
the reference curve (cfr. the classical problem where the approach to the reference 
curve always occurs from below) and also might have a negative initial part. 

In the case in which there is more than one t±, the form of the hump depends on 
the value chosen. In particular, there is no negative part in the hump if we take t\ 
such that £( Ml ~* M2 > (ti) is growing with t. 

Finally, in the glassy phase (T = 0.2 and M = 0.5) the oscillation in S(t) is 
almost completely damped. The hump has a very similar behaviour to the one found 
after a shift in T. Figure IT2l shows the hump for several values of the pairs (M\,ti) 
given in the caption. Also in this glassy case, the height of the hump increases with 
the difference M 2 — Mi and, simultaneously, the position of the maximum has a very 
smooth drift towards smaller values of t — t\. 

Thus, as far as the Kovacs' effect is concerned, we see that the quantum problem 
also shows a non-trivial dependence on the parameter M\ {Ti) and a slow relaxation 
after the perturbation. 
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Figure 11: Hump in the quantum p = 3 model at T = 0.75. The reference curve, 
solid red line, corresponds to M% — 0.5. The modified curves have been obtained 
using, from bottom to top at the maximum, M\ = 0.6, t x — 95 (dashed green); 
Mi = 0.7, t x = 82 (dashed blue); M x = 0.8, t x = 76 (dashed gray); M x = 0.8, t x = 62 
(dashed cyan); M x — 0.8, t x = 46 (dashed magenta); M x = 1, t x = 47 (dashed black) 
and Mi = 1.2, t x = 49 (dashed orange). 
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Figure 12: The hump as a function of the time-difference t — t x in the glassy phase 
of the quantum p = 3 model using, from bottom to top, Mi = 0.6, t x = 539 (dashed 
green), Mi = 0.7, t x = 274 (dashed blue), Mi = 0.8, t x = 204 (dashed magenta), 
Mi = 1, t x = 163 (dashed cyan), Mi = 1.2, t x = 150 (dashed black), and the reference 
curve at constant M2 = 0.5 (solid red). The temperature is T = 0.2 in all cases. 
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5 Conclusions 



We conclude that models with no spatial structure, like the p spin spherical model 
that is intimately related to the schematic mode-coupling theory, can reproduce non- 
trivial memory effects when their non-equilibrium dynamics is studied. Similarly to 
what observed when reproducing the hole-burning protocol |18j . we found here that 
the Kovacs' memory effect is captured by this model. In this sense, the Kovacs' 
experiment is not able to prove the existence of a growing correlation length in 
glassy systems. Assuming that a length scale exists one could, however, compare 
the outcome of this and other experiments to what can be derived from a domain- 
growth like picture for glassy dynamics. 

The reason for having these non-trivial long-memory effects in these fully - con- 
nected spin models is that close to their dynamic critical temperature (and below it) 
their response function after strong perturbation has been applied is not given by a 
simple exponential relaxation. The slow decay of the reponse implies that the effect 
of non-linear perturbations takes very long to disappear. This is encoded in the 
Schwinger-Dyson equations which describe the dynamical evolution of the system. 
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